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(57) Abstract 

The present invention relates to the estimating of CIR (channel impulse response) and SINR (signal-to-interference-plus-noise ratio) 
in CDMA receivers. Particularly, the present invention is directed to a method of determining the channel impulse response (CIR) of a 
communication system, such as the CIR of radio channels of a digital mobile radio network (GSM network). In particular, the present 
invention relates to determining the CIR based on the reception of a known training sequence. The present invention also has application 
in interference cancellation and use in CDMA receivers. The present invention stems from the realisation that the problems associated with 
the prior art can be alleviated by initially transforming the matrix obtained by virtue of network analysis into a 'diagonal matrix* and then 
inverting the resultant diagonal matrix. A further invention is to multiply the reciprocal of the eigenvalue in the inverse by a decreasing 
number until the sidelobes are removed or so small that they have little effect. 
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1 

CHANNEL IMPULSE RESPONSE ESTIMATION USING SINGULAR VALUE DECOMPOSITION 
FIELD OF INVENTION 

The present invention relates to the estimating of CIR (channel impulse 
response) and SINR (signai-to-interference-plus-noise ratio) in CDMA receivers. 
5 Particularly, the present invention is directed to a method of determining the 
channel impulse response (CIR) of a communication system, such as the CIR of 
radio channels of a digital mobile radio network (GSM network). In particular, 
the present invention relates to determining the CIR based on the reception of a 
known training sequence. The present invention also has application in 
10 interference cancellation and use in CDMA receivers. 
BACKGROUND 

In order to determine CIR, a portion of the transmit signal must be known. 
For a GSM network. Synchronisation bursts (SB) are a useful portion of the 
signal. The SB are transmitted on at least one channel from every base station, 
1 5 and they are transmitted in a regular pattern. Decoding of the GSM protocols is 
not necessary. Both the data in the SB and the pattern with which they occur is 
fixed and substantially identical for all base stations. 

The advantage of using SB for determining the CIR is that they represent 
a relatively long, noise-like predetermined transmitted signal. Typically, 64 bits 
20 are transmitted over a period of 237^sec. The determination of the CIR thus 
requires sufficient synchronisation to the bursts in the received signal to enable 
the extraction of the SB which is then processed to determine the CIR. 

The processing is done by using estimation techniques. In order to 
estimate the CIR, a known training sequence is transmitted S tx (t), and this is 
25 corrupted by a communications channel producing the received signal S n (t). 

The problem in estimating the CIR is to determine tap-weights {a} of a FIR filter 
(which is used to approximate the CIR), so that the known Si x (t) after passing 
through the filter is as close as possible to the received signal Srx(t). 

Given the that the transmitted Signal S tx (t) and the received signal S^t) 
30 is known, the CIR can be estimated from: 
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1 . known samples of the transmitted signal T h = S w (t 0 + kt), - N c < k < 
N (note that N + N c samples represent the whole training sequence, the 
numbering has been arranged to start from the known clean samples of the 
received signal which occur N e samples after the start of the training sequence.). 

5 and 

2. measured samples of the received signal Rk = Srx (to + kx), 0 < k < 
N (the first N c samples of the training sequence are assumed corrupted and are 
ignored). 

The tap weights {or} are determined by simple correlation as 



10 



-^Tj-kRi 0<k£N c (2.1) 



j-0 



This algorithm relies on the noise-like properties of the transmitted signal 
whereby its autocorrelation function should have low time sidelobes. The 
sidelobe performance of this algorithm, however, has been found to be limited 

1 5 due to the correlation properties of the -clean" part of the training sequence and 
the fact that only partial correlations are performed for later weights. 

US 5,473,632 estimates CIR by the use of a "least square approach", 
where the Complex Pulse Response (CIR) is predicted by minimising the square 
of the error between what the channel actually outputs and what the estimated 

20 predicted model outputs. This approach gives rise to the equation at Column 6. 
line 35 of US 5, 473,632 which is known from various text books. 

There is considered to be a problem with the US 5,473,632 approach. 
Basically, US 5,473.632 produces a matrix. The problem really relates to what 
is done to the matrix to calculate the CIR. 

25 In US 5.473,632. referring to column 5 lines 27 to 61 and claim 2 thereof, 

a "factor a" is added to the matrix before the matrix is inverted, and then after the 
"factor a" is added, the matrix is inverted. In practice it has been found that the 
inversion of the matrix causes problems in calculating the CIR. 
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SUMMARY OF INVENTION 

The present invention seeks to alleviate the problems experienced in 
determining CIR and SJNR with prior art techniques. 

In essence, the present invention stems from the realisation that the 
5 problems associated with the prior art can be alleviated by initially transforming 
the matrix obtained by virtue of network analysis into a 'diagonal matrix 1 and 

then inverting the resultant diagonal matrix. Thus there is no lactor a* added. 

An advantage of using a diagonal matrix is that it is possible to remove 
or ignore some eigenvalues (such as some - of the small problematic 
10 eigenvalues) and in this way avoid the magnified errors and noise of US 
5,473,632 when the matrix is inverted. 

Another advantage of using a diagonal matrix is attributable to the nature 

of a diagonal matrix and that is, the matrix can be inverted by inverting only the 
diagonal elements. This exposes the difficulty caused by small eigenvalues - they 
-| 5 lead to very large elements in the inverse. 

However, by ignoring or removing some eigenvalues a further problem 
may develop as shown in figure 1 , that is the production of large sidelobes. 

Thus, the present application discloses a still further invention which is 
directed to alleviating the problem of the production of large side lobes. 
20 This further invention is considered an improvement on the invention noted 
above and is based on, not discarding small eigenvalues (which is the 
approach noted above) but instead to multiply the reciprocal of the eigenvalue 

in the inverse matrix by a decreasing number until the sidelobes are removed or 
so small that they have little effect. 
25 DETAILED DESCRIPTION AND BEST MODE 

A preferred embodiment of the inventions disclosed above will now be 
described with reference to the accompanying drawings, in which: 

Figure 1 illustrates CIR estimated using the algorithm disclosed in US 
5,473,632 using e • 200, S/N = OdB, 
30 Figure 2 illustrates CIR estimated using the algorithm disclosed in US 

fc 5,473,632 using e = 400, no noise, 

Figure 3 illustrates CIR estimated using SVD and including 10 largest 
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eigenvalues, 

Figure 4 illustrates CIR estimated using SVD and including 16 largest 
eigenvalues! 

Figure 5 illustrates CIR estimated using SVD and including 25 largest 
5 eigenvalues. 

Figure 6 illustrates CIR estimated using SVD and including 36 largest 
eigenvalues, 

Figure 7 illustrates OR estimated using SVD and using ramping function 
with ^=4, 02=25 and no added noise, 
10 Figure 8 illustrates CIR estimated using SVD and using ramping function 

with n 1 =4, n a =25 S/N - QdB. 

Figure 9 illustrates CIR estimated using SVD and using ramping function 
with ^=4, n 2 =20 no noise, 

Figure 10 illustrates CIR estimated using SVD and using ramping 
1 5 function with n n =4, n 2 ~20 S/N = OdB, 

Figure 11 illustrates CIR estimated using SVD and using ramping 
function with n<i=1, n 2 =25 no noise, 

Figure 12 illustrates CIR estimated using SVD and using ramping 
function with ni«1 , n 2 =25 S/N = OdB, 
20 Figure 13 illustrates CIR estimated using SVD and using ramping 

function with n 1= =1 . n^=64 no noise (Note dB scale to -6QdB), 

Figure 14 illustrates CIR estimated using SVD and using ramping 
function with n n =1, n 2 =64 S/N - 40dB, 

Figure 15 illustrates CIR estimated using SVD and using ramping 
25 function with n 1= 1 , n 2 =25 S/N = 40dB, 

Figure 16 illustrates CIR estimated using SVD and using ramping 
function with m=1 . n 2 =25 Path 1 at 6.9>is OdB, Path 2 at 1 8.5^s -20dB, Path 3 at 
29.5jis -10dB no noise, 

Figure 17 illustrates CIR estimated using SVD and using ramping 
30 function with n n «1, n^S Path 1 at 6.9^s -10dB. Path 2 at 18.5^s OdB, Path 3 at 
29.5*is OdB, Path 4 at 40.4ps -6dB S/N - 20dB. and 

Figure 18 illustrates CIR estimated using SVD and using ramping 



WO 98/38772 PCT/FI98/00159 



function with ni=1 , n2=25 Path 1 at 6.9jis QdB, Path 2 at 12.2}is OcJB no noise. 

One method of estimating the tap weights is to determine the weights 
which predict a received signal most closely matching the measured signal. 

Given the transmitted Signal Su(t) with known samples: 

5 T k -Sa(to+kt) ( -M<k<N 

(M is the assumed extend of the CIR - the number of "dirty" samples) 
and a received Signal Snc(t) with measured samples: 
R^S^to+ktJ.O^^N 

The desired FIR tap weights (o^) are estimated by requiring that the 

10 cumulative square error between the measured received signal and the 
estimate of the received signal (i.e. the known transmitted signal passed through 
the FIR approximating the channel) is minimised, 
i.e. minimise 



ivm 
k=o 



M-1 
j«o 



(2.2) 

15 this may be reduced to the following system of simultaneous equations 
Z( a iETn.jT n ,)=i:R k T M 0<i<M (2.3) 

j.O \ k-0 / k»0 

or using matrix notation 

Ac = TR (2.4) 
where the matrix A is the MxM correlation matrix given by 

ai,i=2T n .^i (2-5) 

20 ^ 

and T is the MxN matrix of conjugated delayed signals 

tij-Ty Q£i£M 0^j<N (2.6) 

and vector c is the required CIR (cxk) 

Note that simple correlation corresponds to 
25 c = TR (2.7) 

The matrix equation (2.4) has the formal solution 



WO 98/38772 PCT/FI98/00159 

6 

c = (A-nT)R (2.8) 
showing that the CIH. may be determined from the samples of the received 
signal using a simple matrix multiplication, 
f p a *t Squares Fstfmation using US 5.473.632 
5 Equation (2.8) equates to the equation disclosed in US 5.473,632 at cot 

6, line 35, 

The complication with the use of this equation arises when the equations 
are solved directly in that they are generally poorly conditioned. In fact, it has 
been found that for the GSM training sequence in the SB the eigenvalues of the 
10 matrix vary over a range exceeding 1010. This is considered to make the 
inverse of the matrix extremely susceptible to roundoff error and the result 
extremely susceptible to noise on the measured signal. 

It has been found that by applying conventional techniques to invert this 
matrix or to solve the equations, the results are not considered useful, even 
15 when using IEEE double precision arithmetic. 

US 5,473,632 discloses an algorithm which stabilises the inversion of the 
matrix and produces results reasonably resilient to noise. US 5,473,632 
estimates the CIR by 

c-((A + £l)-iT)R (2.9) 
20 where I is the identity matrix and e is a constant (they call a noise term). This 
equation equates to the equation disclosed in US 5,473,632 at col. 6, line 44. 
The results of using this algorithm follow: 

The performance of this algorithm critically depends on the choice of the 
variable e. For the matrix A which has diagonal elements with magnitude 192, 
25 the following figures show the performance of US 5,473,632 algorithm using 
difference values of e in determining a CIR consisting of a single path. 

Referring to Figures 1 and 2,for example, it can be seen that there is a 
trade-off between the resolution (width of the main lobe) and the noise 
performance. The sidelobe level of approximately -18dB seems a characteristic 
30 of this algorithm. 

We achieved a best comparable performance by using e= 100, and on 
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this basis comparisons with US 5,473,632 are made. 
Least Squares Estimation using SVD 

In the present invention, we have chosen to apply the use of Singular 
Value Decomposition (SVD) to the poorly conditioned matrix for performing 
5 least squares estimation. In this way, it has been found that we transform the 
matrix that is to be inverted to a diagonal matrix and then the new created 
diagonal matrix is inverted. 

Using the SVD the matrix A (which is positive definite) may be 
represented by 
10 A = USV« (2.10) 

where U and V are orthogonal and S is diagonal. In fact the diagonal elements 
of S are the eigenvalues (all positive) of A and these are arranged in decreasing 
order: 

S 1p i > S 2f 2 ^S ntfl (2.11) 

1 5 If the matrix is singular, some of these are zero. The inverse of A is then 

A-1 =VS-1U< (2.12) 
where S- 1 is formed by inverting the (diagonal) elements of S. The effects of a 
poorly conditioned matrix are clear here, the (almost zero) eigenvalues are 
inverted and produce an extremely large contribution to the inverse. These 
20 amplify the noise. 

Traditional SVD techniques only invert some of the eigenvalues and stop 

at some point. Instead of placing 1/A.t into the inverse of S, the present invention 
uses zero instead. The performance of these techniques along with others are 
evaluated below. 

25 The coefficients in the matrix in equation (2.8) are purely functions of the 

known training sequence, hence they can be computed beforehand. This matrix 
can be inverted once and then the solution at any time requires a matrix 
multiplication. Thus to obtain the CIR estimate from a set of signal samples 
requires M * N complex multiplications and M * (N-1) complex additions, a total 

30 of 6*M-N + 2*M* (NM ) operations. 

The application of SVD to the overdetermined least squares problem 
requires only including the contribution from the N largest eigenvalues in the 
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pseudo-inverse of the matrix, in the case being modelled, the matrix has 
dimension 64, hence N is in the range 1 to 64. The following plots show the 
performance of the SVD algorithm for N- = 10, 16, 25, 36. The effect of 
increasing the number of eigenvalues is to increase the resolution of the 
5 estimator, that is reducing the width of the main fobe. 

Although not shown, simulations were performed which showed that the 
noise performance degrades rapidly in the range of 10 to 20 eigenvalues. 

The present SVD approach should yield good performance in noise as 
the rank of the pseud-inverse is minimised. The most bothersome aspect of the 

10 results is the persistent time sidelobes at approximately -13dB. By analogy with 
Fourier theory, where a truncated time waveform will give frequency sidelobes 
of -13dS f it was discovered that the ramping down of the contribution from 
successive eigenvalues instead of simply cutting off after a certain number 
proved very beneficial in providing control over the time sidelobes. 

1 5 A number of ramping functions may be used such as Bartlett, Blackman- 

Harris, Poisson, Riesz, Hamming, Kaiser-Bessel, Lanczos, Tukey and other well 
known windowing functions. The function chosen in this embodiment has been 
derived from the Hanning window although the raised cosine ramp defined as 
below proved most successful: 




(2.13) 



where w$ weight applied to the contribution from the i th eigenvalue. As an 
example. Figure 7 shows the effect of n 1 =4, n 2 =>25. 

Note that the -l3dB sidelobes have been reduced to approximately- 
22dB (already significantly better than the performance of US 5,473,632). The 
25 width of the main lobe is narrower than the e=1 00 US 5,473,632 algorithm with 
lower sidelobes, The performance with a S/N ration of OdB is shown in figure 8. 
This performance is similar to that achieved with the e^lOO US 5,473,632 

solution, however note that the main lobe is narrower, this will degrade this 
noise performance. 
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Noise performance can be improved at the expense of the lobe width, as 
illustrated in Figures 9 and 10 which have n^, n2=20. 

Considerable trials have taken place to find a satisfactory solution. The 
goals were: 

5 • greater than 20dB dynamic range (preferably 25dB) to give a 

significant edge over the performance of US 5.473,632 
noise performance similar to correlation at S/N - OdB. (This says 
that we don't want the solution to unnecessarily magnify the noise 
contributions). 

10 The results determined from the above uses the weighting factors 

obtained with n^l. n 2 =25 and the performance is shown in figures 11 and 12. 

Performance of the present algorithm on more complicated CIR's will be 
shown shortly. It is interesting to see how much dynamic range the ramping 
algorithm can provide as illustrated in the example shown in Figure 13. 

15 This example is extremely sensitive to noise. Figure 14 shows the 

performance with S/N = 40dB. Note that even this small amount of noise is 
greatly magnified. This same amount of noise produces no perceptible noise on 
the preferred results as determined above and as shown in Figure 15. 

Figures 16 and 17 illustrate the results for some more complex CIR's. 

20 It has also been found that the present embodiment has the ability to 

enable the algorithm to discriminate between peak, and Figure 18 shows close 
peaks split with time separations greater than about 5ms. 

Matrix inversion as described above is also required when estimating the 
channel impulse response vector h in context of interference cancellation 

25 Oh=Y (2.14) 

where <J> is the crosscorrelation matrix between locally generated signals and ¥ 
is the crosscorrelation matrix between locally generated signals and the 
received signal.' The principles and method as described above are equally 
applied in this situation. 

30 Similarly, in interference limited environment CDMA is applied usually 

using optimum combining of signals received from multiple sensors (antenna 
array). Optimum combining maximises the signal-to-interference-plus-noise 
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ratio (SINR). The weight vector for the antenna array to maximise SINR is 
obtained, for example from 

w»aR-i Ud* • (2.15) 
where w is the weight vector for antenna elements, a is a constant, R is the 
5 received interference-plus-noise correlation matrix, and u d * is the conjugate of 
the desired signal vector* We have found an efficient way of getting the weight 
vector by first calculating the CIR. Some optimum combining techniques 
calculate first the CIR (and the noise covariance matrix Q) which is then used in 
multidimensional MLSE to combine signals from the antenna array. Usually the 
1 0 matrix inversion is problematic in mobile radio channels. The embodiment 
as described above can be applied directly to matrix inversion. 
To obtain the weight vectors or employ the inverse of Q to optimally combine the 
signal. 

In summary, whilst alternatives and variants are intended to be 
15 encompassed by the present application, the preferred implementation of the 
inventions disclosed is to generate a CIR using the matrix multiplication as 
specified in equation (2.8) and to then use SVD to determine the inverse of the 
matrix A. The eigenvalues included in the inverse are ramped (weighted) 
according to equation (2.13) with ^=1, n2»25. 



WO 98/38772 



-11- 



PCT/FI98/00159 



CLAIMS: 

1 . * A method of estimating a channel impulse response (CIR) of a 
communication system for use by an estimated predicted model of a 
communication system, the method comprising the steps of: 

predicting the CIR by minimising the square of the error between a 
measured received signal and a predetermined estimate of said received signal, 
said estimate being based on an approximation of a radio channel; 

generating a first matrix form of said CIR; 

transforming said first matrix to a diagonal matrix; and 

inverting said diagonal matrix. 

2. A method as claimed in claim 1 , whereby the diagonal matrix is 
obtained by applying Singular Value Decomposition (SVD) to the first matrix. 

3. A method as claimed in either claim 1 or claim 2 t comprising the 
further step, following the step of inverting, of ramping down the contribution of 
excessive eigenvalues occurring in the inverted diagonal matrix. 

4. A method as claimed in claim 3, whereby a windowing function is 
used to ramp down excessive eigenvalues. 

5. A method as claimed in claim 4 f whereby any one of Bartlett, 
Blackman, Harris, Gaussian, Poisson, Riesz, Hanning (raised cosine), Kaiser- 
Bessel, Lanczos, Tukey or other similar known windowing functions are used to 
ramp down said excessive eigenvalues. 

6. A method as claimed in any one of the preceding claims, when 
used in a GSM network, interference cancellation or CDMA receivers. 
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